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ABSTRACT 

The persistent Be/X-ray binary RX J0440.9+4431 flared in 2010 and 201 1 and has been followed by various X-ray facilities {Swift, 
RXTE, XMM-Newton, and INTEGRAL). We studied the source timing and spectral properties as a function of its X-ray luminosity to 
investigate the transition from normal to flaring activity and the dynamical properties of the system. We have determined the orbital 
period from the long-term Swift/BXT light curve, but our determinations of the spin-period are not precise enough to constrain any 
orbital solution. The source spectrum can always be described by a bulk-motion Comptonization model of black body seed photons 
attenuated by a moderate photoelectric absorption. At the highest luminosity, we measured a curvature of the spectrum, which we 
attribute to a significant contribution of the radiation pressure in the accretion process. This allows us to estimate that the transition 
from a bulk-motion-dominated flow to a radiatively dominated one happens at a luminosity of ~2 x 10 36 ergs _1 . The luminosity 
dependency of the size of the black body emission region is found to be tbb °c L^ 39±0 02 . This suggests that either matter accreting 
onto the neutron star hosted in RX J0440.9+4431 penetrates through closed magnetic field lines at the border of the compact object 
magnetosphere or that the structure of the neutron star magnetic field is more complicated than a simple dipole close to the surface. 
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1. Introduction 

LS V +44 17/RX J0440. 9+4431 was classified as one of the rare 
persistent neutron star Be/X-ray binaries based on the outcomes 
of the first dedicated RXTE monitoring campaign performed by 
Reig & Roche (1999). The source distance was later estimated 
at (~3.3 ± 0.5) kpc by Reig et al. (2005). RXJ0440. 9+4431 
is also reported in the seven-year and nine-year INTEGRAL 
all-sky survey (Krivonos et al. 2010b, 2012). From the former 
to the latter catalog, the average flux of the source increased 
from (0.95 ± 0.15)mCrab to (4.0 ± O.ljmCrab 1 (17-60keV) 
as it underwent several episodes of enhanced X-ray activity in 
2010 and 2011 (see later in this section). The flux estimate 
given in the seven-year catalog is consistent with what has been 
measured previously by Reig & Roche (1999) and can thus be 
considered as the persistent level of the source. A similar flux 
(0.9 ± 0.5)mCrab 2 was also reported in the 54-month Palermo 
BAT Survey in the 15-150keV band (Cusumano et al. 2010). 

The first observed outburst from RXJ0440. 9+4431, which 
we are aware of, was detected on 2010 March 26, and lasted until 
April 15 (Morii et al. 2010; Finger & Camero-Arranz 2010). The 
Swift / BAT light curve (Fig. 1) showed that the source reached a 
peak luminosity of about 200mCrab (15-150keV) six days af- 
ter the onset of the event. It remained in this state for eight days, 
and then turned off to its persistent luminosity level in about 

1 In the energy band 17-60 keV these correspond to 1.35 and 5.7 x 
10 -11 erg cm -2 s _1 , respectively. 

2 In the energy band 15-150 keV this corresponds to 2.0 x 

10 -11 erg cm -2 s -1 . 


ten days. The outburst was also observed by the RXTE All Sky 
Monitor (2-10 keV, ASM). In the ASM light curve, the source 
reached about lOcts/s (i.e., ~130mCrab). A detailed study of 
the first outburst from RXJ0440. 9+4431 was presented by Usui 
et al. (2012). The authors show that the broad-band spectrum of 
the source could be described well by a combination of a black 
body component, a cut-off power law, and an iron line with a 
centroid energy of ~6.4 keV. They also find a pronounced dip in 
the pulse profile interpreted as being due to obscuration of the 
X-ray emission from the neutron star by the accretion stream. 

A second outburst from RXJ0440. 9+4431 was detected by 
INTEGRAL on 2010 September 1 (Fig. 1). The measured av- 
erage and peak fluxes during the event were ~30mCrab and 
(77 ± 17)mCrab (17-60 keV), respectively (Krivonos et al. 
2010a). The SwiftfBAI and RXTE/ASM light curves showed 
that the outburst had a rise and decay time of about five days, 
with a compatible peak luminosity of ~80mCrab. The analy- 
sis of the combined Swift, INTEGRAL, and RXTE data led to 
identifying a possible cyclotron absorption feature at ~30 keV, 
suggesting a surface magnetic field of B - 3 x 10 12 G for the 
accreting neutron star hosted in this system. This magnetic field 
intensity would be compatible with the measured properties of 
the source broad-band noise in the X-ray domain (Tsygankov 
et al. 2012). 

The two intense X-ray emission episodes recorded from 
RXJ0440. 9+4431 showed the typical properties of the so-called 
“Type I” outbursts, usually displayed when the neutron star in a 
Be/X-ray binary approaches the system periastron and interacts 
with the companion equatorial disk. The presence of a neutron 
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Fig. 1. Light curve of RXJ0440. 9+4431 and times of the observations 
reported in this paper. Black points with errors: Swift/BAT daily aver- 
aged light curve. Red vertical bars: times of the RXTE observations. 
Blue horizontal bar: time span of the INTEGRAL data used in this 
paper. Green diamonds: times of the Swift observations used in our 
analysis. The first vertical yellow bar indicates the occurrence of the 
first maximum in the light curve; the others correspond to the times 
of the following maxima predicted by our orbital period determination. 
MJD 55 280, 55 430, and 55 580 correspond to 2010, March 25, 2010, 
August 22, and 2011, January 19, respectively. 


star in RX J0440. 9+443 1 is secured by the detection of pulsa- 
tions in its persistent and outburst emission. The measured pulse 
period increased from ~202s in 1999 (Reig & Roche 1999) to 
~206 s during the latest observations, thus evidencing an aver- 
age spin down trend. The first estimate of the source orbital pe- 
riod was reported by Reig (2011) at ~150d using the orbit-to- 
spin period relation for Be/X-ray binaries. This period would 
be compatible with the recurrence time of the source outbursts 
(150-160 d Morii et al. 2010; Krivonos et al. 2010a; Tsygankov 
et al. 2012). 

A deep (20 ks) XMM-Newton observation of RXJ0440.9+ 
4431 was carried out on 2011 March 18 to study the prop- 
erties of its persistent emission (La Palombara et al. 2012). 
During this observation, the source displayed an average flux 
of 6 x 10 -11 erg cm -2 s _1 (0.3-12 keV) and the corresponding 
spectrum could be well described using a power law plus a 
black body model. The estimated radius for the thermal emis- 
sion component was compatible with the size of a hot-spot on 
the surface of a strongly magnetized neutron star (~10 12 G). 
Pulsations at ~205 s were detected throughout the observation. 
No significant spectral variability with the spin phase could be 
measured, although an indication of hardening was found at the 
pulse minimum. 

In this paper, we report on the last detected outburst from 
the source which occurred on February 2011 (Fig. 1). On this 
occasion, we triggered our pre-approved RXTE target of oppor- 
tunity observations and performed a simultaneous monitoring 
campaign at the softer X-ray energies with Swift/XRT to study 
the evolution of the source spectrum with luminosity and im- 
prove the measurement of the source orbital period. We have 
also re-analyzed all available RXTE, Swift/X RT, XMM-Newton , 
and INTEGRAL (Bradt et al. 1993; Gehrels et al. 2004; Jansen 
et al. 2001; Winkler et al. 2003) observations performed in the 
direction of RX J0440.9+4431 in 2010 and 2011 to compare the 
source emission properties across two orders of magnitude in 
luminosity. The data analysis technique is described in Sect. 2, 
and all the results are reported in Sect. 3. Our discussions and 
conclusions are given in Sects. 4 and 5, respectively. 


2. Data analysis 

We have analyzed RXTE and Swift data using the HEASOFT 
software package (v. 6.12) while all spectral fits have been per- 
formed using XSPEC version 12.6.0u (Arnaud 1996). 

For RXTE, we have followed the suggested analysis threads 3 
and we considered only PCA (2-60 keV, Jahoda et al. 1996) 
data from the PCU2 unit, which was always active. Good time 
intervals (GTI) were created by imposing that at least 30 min 
passed after the satellite exited the SAA, and the elevation on 
the Earth limb was >10°. We additionally imposed an electron 
ratio <0. 15 and modeled the background using the provided es- 
timator v3.8. Events were extracted in the mode GoodXenon 
(time resolution 0.95 ps and 255 energy channels). The instru- 
ment specific response matrices were generated taking into ac- 
count the corrected energy bounds and the gain corrections per- 
formed by the EDS. The average PCA spectrum was obtained 
from the standard2 mode without further rebinning and fit- 
ted above standard channel 4 (~3 keV) up to the higher energy 
at which the source signal was detectable (i.e., about 20keV, 
for short observations or faint source states, higher otherwise). 
Spectra were assumed to have a systematic uncertainty of 0.5%, 
as recommended by the instrument team 4 . 

For Swift/XRT (0.5-1 1 keV), we followed standard analysis 
procedures (Burrows et al. 2005) and considered both window 
timing (WT) and photon counting (PC) data, but limited our 
spectral range to 0.8-10 keV to avoid calibration uncertainties. 
XRT spectra were grouped to have at least 50 counts in each en- 
ergy channel using grppha. We extracted separated spectra for 
data accumulated in each satellite orbit and used the proper ex- 
posure map to account for bad columns, hot pixels, vignetting, 
and PSF losses. 

We extracted the source spectrum from the INTEGRAL 
IBIS/ISGRI (20-100keV, Ubertini et al. 2003) and JEM-X 
(5-30 keV, Lund et al. 2003) data of the 2010 outburst taken 
from 2010, September 1 to September 10 using the Off- 
line Science Software (OS A, v.10) distributed by the ISDC 
(Courvoisier et al. 2003). The JEM-X spectrum is extracted in 
the standard 16 bins defined in OS A, while the IBIS/ISGRI one 
in user defined 62 bins. 

XMM-Newton data have been processed using XMM 
Science Analysis System (SAS) version 12.0.1. We have ex- 
tracted only the EPIC-PN (0.5-10keV, Struder et al. 2001) 
spectrum from a region of 0.6 arcmin centered on the 
source. Background and pile-up corrections were neglected 
(La Palombara et al. 2012). The spectrum was adaptively re- 
binned using the SAS tool specgroup in order to have at least 
25 counts per energy bin and, at the same time, to prevent over- 
sampling of the energy resolution by more than a factor of three. 

3. Results 

3.1. Orbital period 

To search for the source orbital period, we exploited the daily av- 
eraged light curve obtained from the hard X-ray transient mon- 
itor page of Swift/ BAT 5 . These data cover almost continuously 
the time interval MJD 53414-56282 (2005 February 13 - 2012 

3 http : //heasarc . gsfc.nasa.gov/docs/xte/xhp_proc_ 
analysis.html 

4 http : //heasarc . gsfc.nasa.gov/docs/xte/pca/doc/rmf/ 
pcarmf- 11 . 7/ 

5 http : //heasarc .gsfc .nasa. gov/docs/swift/results/ 
transients 
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Fig. 2. Folded Swift / BAT light curve (15-150 keV) of RXJ0440.9+4431 
at the best determined orbital period, (150.0±0.2) days. The inset shows 
the periodogram obtained from the epoch folding technique applied to 
the Swift/BAI light curve. 


December 31), owing to the large field of view of the instru- 
ment and the spacecraft pointing strategy. We first used the 
Lomb-Scargle periodogram method in the range 32-200 days. 
The highest significance period is found at ~ 150 days. We have 
also exploited the epoch folding technique on the same data by 
assuming 32 phase bins and assigning weights of 1/cr 2 to the 
individual bins (here cr is the reported uncertainty on the source 
count rate). The most prominent peak (y 2 = 942) is found also 
in this case at ~ 150 days with evident harmonics at multiple fre- 
quencies (see inset in Fig. 2). We have refined our determination 
of the orbital period by limiting our search with the above meth- 
ods around the most significant peak and using 130 phase bins. 
The best value is found at (150.0 ± 0.2) days. The uncertainty 
at 68% c.l. on this measurement is obtained by bootstrapping 
1000 light curves with count rates randomly distributed within 
the uncertainty of the measurements, and then computing the 
standard deviation of the period determinations. We have folded 
the light curve at this period and obtained the profile shown in 
Fig. 2. This is reminiscent of what is usually observed from 
Be/X-ray binaries, with most of the emission concentrated dur- 
ing the periastron passage (phase 0.6 in the figure). We have fit 
the folded light curve with a constant plus a Gaussian to estimate 
the time of the light curve maximum. The Gaussian is centered at 
orbital phase 0.641 ± 0.001. which can be converted into the time 
of occurrence of the first flare by taking into account the best or- 
bital period: it corresponds to MJD 55 293.2 ± 0.3, i.e., the time 
of the brightest flare observed so far from RX J0440. 9+443 1 . We 
have plot this time interval in Fig. 1 , left section, as a vertical 
shaded region. We have then propagated the uncertainty on the 
orbital period to the following two flares observed in our cam- 
paign and plot them as shaded regions in the other two sections 
of the light curve. 

Interestingly, when limiting ourselves to the time range be- 
fore MJD 55 280 (2010 March 25), the orbital modulation was 
not detectable and no prominent peaks are clearly visible over 
the red noise. This indicates that the X-ray active phase of the 
source began only in April 2010 when it was detected in out- 
burst for the first time, since the availability of modern X-ray 
monitors. 


Table 1. Measurements of the spin period of RX J0440.9+4431 ob- 
tained from all the RXTE observations analyzed in the present work. 


Time 

MJD 

A T a 
ks 

P ■ 

1 spin 

S 

1 cr error 
s 

55 292.35 

l.n 

205.10 

0.14 

55 298.23 

1.24 

205.03 

0.13 

55 301.23 

3.51 

204.84 

0.13 

55 447.02 

3.60 

204.89 

0.17 

55 447.95 

2.61 

204.24 

0.06 

55 449.06 

3.32 

205.51 

0.07 

55 450.24 

2.60 

204.60 

0.11 

55 451.75 

3.67 

205.38 

1.19 

55 452.99 

3.66 

205.26 

1.36 

55 594.16 

3.67 

204.65 

0.29 

55 595.65 

3.66 

203.97 

0.26 

55 596.50 

2.79 

205.42 

0.19 

55 597.55 

3.07 

204.57 

0.25 

55 598.45 

2.95 

205.00 

0.18 

55 599.50 

2.81 

205.99 

0.20 

55 600.55 

3.35 

205.78 

0.25 

55 601.39 

3.59 

205.86 

0.18 

55 602.37 

3.56 

206.29 

0.18 

55 603.49 

2.83 

206.21 

0.24 

55 604.34 

0.97 

205.25 

0.86 

55 605.45 

2.45 

206.08 

0.45 

55 606.49 

1.48 

206.01 

0.47 


Notes. w AT = (r end - T start ) /2. 


3.2. Timing analysis and spin period measurements 

To search for the ~200 s source spin period and its possible vari- 
ations, we performed a timing analysis of all RXTE data (event 
mode GoodXenon). Swift pointings were on average too short 
to perform a similar analysis. We first extracted the background 
corrected PCA light curves in the 3-35 keV energy range (time 
bin 2 s), and barycenter-corrected them using the known opti- 
cal position (Reig et al. 2005). For each observation we removed 
the variability on time scales longer than the expected pulsed 
signal 6 and determined the pulse period with the epoch folding 
technique. 

We obtained the best determined period for all the avail- 
able RXTE observations using an optimal number of phase bins 
(from 16 to 64 depending on the S/N), and estimated the un- 
certainties by running 100 simulations of the input light curves 
with rates estimated from Poisson distributions centered on the 
real measured values. The results are reported in Table 1. Some 
variability in the measured values is evident from the table, but 
the limited number and accuracy of our determinations are not 
sufficient to disentangle unambiguously the effects of the orbital 
motion and/or possible accretion torques. 

3.3. Broad-band spectral analysis 

The availability of quasi-simultaneous INTEGRAL, RXTE, and 
Swift data permitted to investigate the broad-band spectral prop- 
erties of RXJ0440. 9+4431 (note that part of the data were al- 
ready presented: Tsygankov et al. 2012; Usui et al. 2012). The 
broad-band quasi-simultaneous data set is shown in Table 2 with 

6 To achieve this, we computed the sliding average of the light curve 
on intervals of 400 s every 200 s and then interpolated these values at 
the central time of each time bin to subtract them from the light curve. 
The resulting time series was then used as input to the epoch folding 
algorithm. 
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Table 2. Log of observations of the 2010-2011 outbursts used in this paper. 


ID" 

OBS ID 6 

Exp. 

Start 

Stop 

logiolFhix)" 

Flux" 

U d 



ks 

MJD 

MJD 


xlO -9 

xlO 35 


95418-01-01-00 

3.1 

55 292.3306 

55 292.3778 

-8.2412 ± 0.0007 

5.738 

74.766 


95418-01-02-00 

3.3 

55 298.2106 

55 298.2554 

-8.4228 ± 0.0008 

3.778 

49.222 


00031690001 

4.8 

55 299.8793 

55 300.0424 

- 

- 

- 


95418-01-02-01 

3.2 

55 301.2178 

55 301.2611 

-8.7541 ±0.0013 

1.762 

22.954 


95418-01-03-01 

1.0 

55447.0197 

55 447.0311 

-8.738 ± 0.002 

1.827 

23.804 


IBIS 

56 

55441.2576 

55 448.1785 

- 

- 

- 


JEM-X 

64 

55 441.2576 

55 448.1785 

-8.78 ± 0.01 

1.66 

21.6 


00031690003 

5.7 

55443.0294 

55 443.2021 

- 

- 

- 


95418-01-03-00 

2.8 

55447.9302 

55 447.9728 

-8.894 ± 0.002 

1.275 

16.615 


95418-01-04-00 

2.5 

55449.0456 

55 449.0763 

-9.021 +0.003 

0.953 

12.417 


95418-01-04-01 

1.5 

55450.2285 

55 450.2463 

-9.316 + 0.006 

0.483 

6.296 


95418-01-04-02 

1.1 

55451.7385 

55 451.7563 

-10.37 ± 0.06 

0.042 

0.55 2 


95418-01-04-03 

1.0 

55452.9830 

55 452.9963 

-10.47 ± 0.07 

0.034 

0.444 


96364-01-01-00 

1.3 

55 594.1491 

55 594.1680 

-8.937 ± 0.003 

1.155 

15.051 

1 

96364-01-01-01 

1.5 

55 595.6400 

55 595.6574 

-9.120 + 0.004 

0.759 

9.888 

00031690005" 

0.99 

55 595.6759 

55 595.6875 

- 

- 

- 


96364-01-02-00 

3.2 

55 596.4802 

55 596.5271 

-9.526 ± 0.006 

0.298 

3.882 

z 

00031690005" 

1.46 

55 596.5314 

55 596.5484 

- 

- 

- 


00031690006" 

0.86 

55 597.4755 

55 597.4855 

- 

- 

- 

3 

96364-01-03-00 

3.0 

55 597.4921 

55 597.5728 

-10.33 ± 0.03 

0.046 

0.604 


00031690006" 

1.08 

55 597.5394 

55 597.5520 

- 

- 

- 


96364-01-04-00 

3.2 

55 598.4385 

55 598.4848 

-10.38 ± 0.02 

0.042 

0.547 


96364-01-05-00 

2.9 

55 599.4872 

55 599.5304 

-10.41 + 0.04 

0.039 

0.504 


96364-01-06-00 

2.8 

55 600.5345 

55 600.5763 

-9.362 ± 0.004 

0.434 

5.659 


00031690009 

1.62 

55 600.5450 

55 600.5638 

- 

- 

- 


96364-01-07-00 

3.3 

55 601.3754 

55 601.4217 

-10.40 ± 0.05 

0.040 

0.515 


96364-01-08-00 

3.3 

55 602.3547 

55 602.4000 

-10.40 ± 0.04 

0.040 

0.517 


96364-01-09-00 

3.0 

55 603.4719 

55 603.5132 

-10.54 + 0.03 

0.029 

0.376 


96364-01-10-00 

3.2 

55 604.3171 

55 604.3582 

-10.63 ± 0.06 

0.023 

0.304 


96364-01-11-00 

3.1 

55 605.4297 

55 605.4710 

-10.63 ± 0.03 

0.023 

0.302 


Notes. The three vertical blocks separated by a horizontal line refer from top to bottom to April 2010, September 2010, and February 2011, 
respectively. <a) ID is the conventional name used in our broad-band spectral analysis exploiting quasi-simultaneous observations. ® OBS IDs as 
9xxxx-01-xx-xx refer to RXTE, 0003169000x to Swift, IBIS and JEM-X to the INTEGRAL monitoring. (c) 2-30 keV absorbed flux is in units of 
erg cm -2 s -1 , not reported for instruments operating in different bands. w 2-30 keV luminosity in erg/s assuming a distance of 3.3 kpc. (e) Different 
snapshots with the same OBS ID. 


a conventional ID of the corresponding multi-instrument obser- 
vations: A, S, 1, 2, 3, 4. For the INTEGRAL data reduction 
we took advantage of the most updated IBIS/ISGRI calibrations 
included in the OS A v.10 7 . 

To compare our findings with those reported previously by 
Tsygankov et al. (2012), and test the effect of the improved 
INTEGRAL calibrations, we used the same data set as in 
the above paper (INTEGRAL data from MJD 55 440-55 449, 
RXTE/PCA data from OBSID 95418-01-03-00 and SwiftfKRT 
pointing OBSID 00031690003, i.e., observation S in Table 2). 
We restricted the selection of pointings to a maximum off-axis 
angle of 7° to minimize the systematic uncertainties on the IBIS 
instrumental response. The total effective exposure time of the 
final IBIS/ISGRI spectrum was 56 ks. This corresponds to an ef- 
fective exposure time for JEM-X2 of 64 ks (owing to the lower 
dead- time of the detector). 

A spectral model comprising a cut-off power law and 
a black body component could fit the data reasonably 
well (^ 2 /d.o.f. = 1.09/558). Normalization constants (fixed 
to 1 for Swift/X RT) were also included in the fit to ac- 
count for source variability and instrument inter-calibrations. 

7 http : //www. isdc .unige . ch/newsletter . cgi?n=n23 


Photoelectric absorption is accounted for using the model phabs 
with abundances from Anders & Grevesse (1989). Some resid- 
uals are present at >30keV, as reported by Tsygankov et al. 
(2012). However, we found that an additional Lorentzian absorp- 
tion feature with cr = 6 keV at ~30 keV, as proposed by these au- 
thors, did not significantly improve the fit (^ 2 /d.o.f. = 1.05/556, 
33% chance improvement probability). 

A better description of the spectrum could be obtained 
by using the general Comptonization model BMC (Titarchuk 
et al. 1997) with an exponential cutoff at high energy (BMC x 
exp(-£ , /£ , Fo id), ^ 2 /d.o.f. = 0.92/559). Because of the general 
shape of the convolutional kernel, this model can be applied to 
Comptonization spectra of both optically thin and optically thick 
regimes, and considering the cases in which thermal or thermal 
plus bulk Comptonization dominate (Farinelli et al. 2013). The 
emerging spectrum of the BMC model can be expressed as 

£ ( 1 ) 

1 + A 

where G(E, Eo) is the Green’s function of the Componization 
energy operator, and S (E) is the seed photon black body spec- 
trum. The quantity A /(A + 1) empirically provides the fraction 
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Table 3. Best fit parameters for the absorbed BMC with exponential cut- 
off model for the broad band spectral fits of the 2010 outburst (uncer- 
tainties are at 90% c.l.). 



April 2010“ 

September 2010* 

A h [ 10 22 cm -2 ] 

0.31+0.02 

0 27 +0 03 
- 0.02 

FfoU [keV] 

29/ 4 

53 +16 
J -10 

kT BB [keV] 

1.38 + 0.02 

1 "57+0.03 
-0.04 

a 

0.50 ± 0.07 

0.78 + 0.12 

l°gl0 A 

>2.15 

0.833 

Abmc [^] c 

0 030 + ° 003 
U,UJU _0. 001 

0.019 ± 0.002 

r BB ‘' fkm] 

2.64 + 0.12 

2.13 + 0.16 

^Gauss [keV] 

zr 0-3+0.04 
-0.02 


^ Gauss 

(6.6+ 1.3) x 10“ 4 


CpcA 

1.55 + 0.02 

0.632 ± 0.008 

ClSGRI 


0.73 ± 0.03 

CjEM-X2 


0.80 ± 0.02 

FluX 2 -30keV e 

3.78 ± 0.01 

1.272 + 0.009 

FluXo.l- 200 keV^ 

3.4 

2.9 

T?ed/ d -°- f - 

1.24/571 

0.93/559 


Notes. M Observation A in Table 2. ® Observation S in Table 2. (c) L39 
is the X-ray luminosity in units of 10 39 erg/s, Dio the source distance in 
units of 10 kph . td> Quantity derived from N BMC and kT BB . (e) Flux of the 
RXTE/PCA observation in units of 10 -9 erg cm -2 s _1 . kf> Estimated de- 
absorbed flux in units of 10 -9 erg cm -2 s _1 , normalized to the SwiftPCKT 
observations. 


of seed photons which are Compton scattered; it is related to 
the geometrical configuration of the system and to the seed pho- 
ton spatial distribution in the Comptonization region. Actually, 
log 10 (A) is the model parameter returned by XSPEC. The Green’s 
function G{E , £ 0 ) is a broken powerlaw oc E a+3 for E < Eq and 
oc E~ a for E > E 0 . The value of the powerlaw energy index a 
is related to the efficiency of the Comptonization process. From 
the normalization of the model, it is possible to infer the effec- 
tive radius of the region in which the seed black body radiation 
is generated by using the relation (Vbmc = T39 /D 2 0 (here L39 is 
the luminosity in units of 10 39 erg/s and Dio the distance in units 
of lOkpc): 


>BB,km 


7.75 x 10 3 


■Dio^bmc 

C(kT)l y 


( 2 ) 


Here C = 1 for a spherical geometry and C = 0.25 for a circular 
slab. Assuming for RX J0440. 9+443 1 a distance of 3.3 kpc, the 

relation above can be written as reB.km = ■\j 844 Nbmc/ C( kT)^,y. 

The best fit parameters obtained using the absorbed BMC 
with exponential cut-off model are reported in Table 3. We found 
that this description of the source continuum emission permit- 
ted to constrain the optical depth of the previously suggested 
Lorentzian shaped absorption feature to <0. 1 (see the residuals 
from the best fit in in Fig. 3). 

To test further the applicability of the spectral model in- 
troduced above to RXJ0440. 9+4431, we carried out a second 
fit to the quasi-simultaneous Swift /X RT and RXTE/PCA obser- 
vations performed in April 2010 (observation A in Table 2). 
At this epoch, the source was brighter (F2-30kev — 3.8 x 
10 -9 erg s _1 cm -2 ) than in observation S, PCA data could be used 



Energy [keV] 

Fig. 3. Broad-band spectrum of the quasi-simultaneous 
INTEGRAL/ISGR1 (green), INTEGRAL//EM-X2 (blue), RXTE/ 
PCA (red), and Sw//f/XRT (black) observations performed in 
September 2010. The best fit is obtained with the absorbed BMC with 
exponential cut-off model. The residuals from this fit are also shown in 
the bottom panel. 


from 3 to 45 keV and an iron line was needed to fit the data 
(we modeled this as a thin additional Gaussian component cen- 
tered at ~6.4keV; Usui et al. 2012). No evidence was found for 
the Lorentzian shaped absorption feature (width 6 keV, centroid 
energy <42keV) suggested by Tsygankov et al. (2012). We ob- 
tained a 90% c.l. upper limit of 0.03 on its optical depth. We also 
note that Usui et al. (2012) did not find evidence of any absorp- 
tion feature in the combined RXTE/PCA and HEXTE spectrum 
of the source. We conclude that the absorption feature reported 
by Tsygankov et al. (2012) was probably due to systematic bi- 
ases in the joint analysis of data from several instruments and/or 
a inadequate modeling of the continuum emission. 

3.4. Spectral variability 

To study possible changes in the spectral parameters of the 
source at different luminosity states, we fit here all the available 
RXTE/PCA data, which provided an unbiased monitoring of the 
source across all different states, by assuming the continuum 
model for the source broad-band emission discussed in Sect 3.3. 
A log of the observations is given in Table 2, where we report 
also the source flux and luminosity. The lower energy thresh- 
old of the RXTE/PCA is 3 keV and it does not permit a precise 
measurement of the absorption column density in the direction 
of RX J0440. 9+443 1 . Thus, we have fixed it to the weighted 
average value measured from all the available Swift/XRT and 
RXTE quasi-simultaneous observations characterised by high 
enough signal-to-noise ratio (S/N) to significantly constrain our 
parameter of interest (the observations with a conventional ID 
in Table 2). We find Ah = (2.8 ± 0.1) x 10 21 cm -2 . In several 
fits the value of log 10 (A) in the BMC model could only be poorly 
constrained. Where possible, we estimated a lower limit at 90% 
c.l. by using the XSPEC command steppar. For the lower statis- 
tic spectra we fixed it to log 10 (A) = 1.5. During the observations 
performed when the source reached the highest flux (observation 
IDs 95418-01-01-00 and 95418-01-02-00, April 2010) an expo- 
nential cut-off at high energy (i.e., exp (-D/Dfow)) was needed 
to properly fit the data: we have found DfoM = (30 ± 6) keV and 
(33 ± 8) keV, respectively. For all the other observations, we have 
verified that the lower limit on the cut-off energy is outside the 
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- 10.5 - 10.0 - 9.5 - 9.0 - 8.5 

log 10 Flux (2-30 keV) [erg s _1 cm~ 2 ] 


Fig. 4. Results of the spectral fits to the RXTE/PCA observations (filled circles) performed during the outbursts of RX J0440.9+4431 occurred 
in 2010 and 2011 (Table 2). The best fit spectral parameters obtained from the XMM-Newton observation (2011 March 18) are also included 
and represented with a star (uncertainties are at 68% c.l. and lower limits at 90% c.l.).The BMC model was used for all fits, but during the two 
observations performed when the source reached the highest flux, an exponential cut-off at high energy was needed to properly fit the data (see 
Sect. 3.4). The black body radius r BB is obtained from Eq. (2) setting C = 1. The red solid lines in the right panels represents the best fit relations 
of Eq. (3), (4) obtained from the PCA data only. 


energy range covered by our data (~50 keV), consistently with 
the value derived using INTEGRAL/ISGRI data (see Sect. 3.3). 
The results obtained from all the fits are reported in Fig. 4. We in- 
clude in the figure also the EPIC-PN data from the XMM-Newton 
observation that caught RXJ0440. 9+4431 in its persistent low 
state (La Palombara et al. 2012). Note that the BMC model also 
provided a fully acceptable fit to these data (^ 2 /d.o.f. = 1.1/183) 
and yielded log 10 A > 1.2, VT = (1.02 ± 0.03) keV, a = 
0.81 ± 0.08, /"bb = (0.73 ± 0.04) km, and a column density 
Nu = (5.4 ± 0.2) x 10 21 cm -2 (uncertainties at 90% c.l.). The ab- 
sorption column density measured with XMM-Newton/EPIC-PN 
was higher than that measured with Swift/X RT. We verified that 
fitting the PCA data with this higher value of would not affect 
the results discussed below. 

From the plots in Fig. 4, we note that the parameter 
log 10 (A) could be only poorly constrained and a shows a not 
very pronounced downwards trend for increasing source flux. 
Conversely, both the effective radius and temperature of the ther- 
mal emission component increase with the source flux. For the 
effective radius, we find (^ 2 d /d.o.f. = 0.6/19): 

log 10 r BB = (0. 1 82 ± 0.008) + (0.39 ± 0.02) (log 10 F_ 9 ) , (3) 

where r B B is the effective radius of the black body seed pho- 
tons expressed in km and F - 9 is the 2-30 keV flux expressed in 
units of 10 _ 9 ergcm _ 2 s _1 . This relation is equivalent to /- B b 
^•a39±o.02 or K r 2.56±o.i3 similarly for the black body tem- 
perature, we find (y 2 ed /d.o.f.=0.6/19): 

log 10 kT ke y = (0.128 ± 0.004) + (0.070 ± 0.008) (log 10 F_ 9 ) ,(4) 

where kT^y is the black body temperature expressed in keV. 
In Eqs. (3), (4) uncertainties are given at lcr c.l. To prove the 
variations in the source spectral properties in the different obser- 
vations in a model-independent way, we also report in Fig. 5 the 


<3> 



0.5 0.6 0.7 0.8 0.9 '.0 1.1 

(6-10 keV)/(3-6 keV) 


Fig. 5. Hardness-intensity diagram of all the available PCA observation 
of RXJ0440. 9+4431. Different symbols refer to different outbursts: di- 
amonds are for the outburst in April 2010 (MJD55 292-55 301), stars 
for the outburst of September 2010 (MJD 55 447-55 450), and black 
dots for the outburst of February 201 1 (MJD 55 594-55 606). 


hardness-intensity diagram of the RXTE/PCA data (as shown by, 
e.g., Reig & Nespoli 2013). From this figure we note that dur- 
ing all outbursts the source X-ray emission becomes harder at 
higher fluxes and there is an evident change of the overall slope 
at PCU2 count rates >100 cts/s (the persistent emission level of 
RXJ0440. 9+4431 corresponds to a PCU2 count rate of about 
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1-3 cts/s) We verified that different choices of the energy bands 
would not significantly affect this result. 

An interpretation of the spectral changes in RXJ0440.9+ 
4431 is provided in Sect. 4.3. 

4. Discussion 

4.1. The long-term X-ray variability and the orbital period 

The Swi/t/BAT data allowed us to measure with a fairly good 
accuracy the orbital period of the source owing to its enhanced 
X-ray activity in 2010 and 2011. Before this period, the avail- 
able data revealed that RX J0440. 9+443 1 was a persistent X-ray 
emitter at a level of about 1 mCrab in hard X-rays (15-150keV) 
and no significant orbital modulation could be detected from 
its light curve. Be/X-ray systems are rarely persistent X-ray 
sources, but many of them showed in the past only sporadic 
episodes of enhanced X-ray activity. The latter is thought to be 
originated by variations in the extension of the Be star equatorial 
disk (see, e.g., Reig 201 1). The existence of persistent Be sys- 
tems supports the idea that the fainter long-term X-ray activity 
of these sources is most likely related to the accretion of the NS 
from the tenuous wind of the companion, while the outbursts 
occur as a consequence of the interaction at the system peri- 
astron between the compact object and the Be equatorial disk 
(see also, Rothschild et al. 2012, for the quiescent emission of 
A 0535+26). From the shaded regions of Fig. 1, we note that 
the maxima of the second and third flares do not correspond 
exactly to a multiple of our best candidate orbital period: one 
lies slightly before, the other after. In Be/X-ray binary systems, 
some variability of the orbital phase occurrence of outburst is ex- 
pected, but in absence of orbital ephemeris based on pulse period 
modulation, we cannot speculate further on this subject. 

4.2. The broad-band spectral model 

Previous works on persistent Be pulsars (Reig & Roche 1999) 
showed that their spectra can usually be well described by a 
model comprising a power law component and a black body with 
an emitting radius of a few hundreds meters (La Palombara & 
Mereghetti 2006, 2007; La Palombara et al. 2009, 2012). These 
sources are characterized by a spin period of the order of hun- 
dreds seconds, orbital period of hundreds of days and a persistent 
luminosity of 10 14-35 erg s -1 . The direct black body emission is 
commonly interpreted as being due to the thermal mound at the 
base of the neutron star accretion column (with a size propor- 
tional to the width of the accretion stream), while the power law 
is the signature of Comptonization of the thermal photons by the 
in falling plasma. Following this interpretation and using data 
from several X-ray observatories, we have verified that also the 
broad-band spectrum of RX J0440. 9+443 1 could be reasonably 
well described by this phenomenological model, although our 
wide energy range requirers the introduction of an exponential 
roll-over at high energies. 

A more insightful description of the parameters regulating 
the X-ray emission from the source can be obtained by using 
the BMC model. This permits to describe in a self-consistent 
way the process in which part of the thermal seed photons are 
Comptonized and up-scattered to higher energies, giving rise to 
a power-law shaped emission with part of the photons escaping 
the Comptonizzation to reach directly the observer. The model is 
a general solution of the process of Comptonization of a spher- 
ically symmetric in-falling plasma but has the major limitation 
that the high energy roll-over, due to either the thermal motion 


of electrons or the deviation from a free-fall velocity profile of 
the moving plasma, is not taken into account. This can be empir- 
ically introduced by means of an additional exponential cut-off 
(exp(-£/£ F oid), Farinelli et al. 2013), still preserving the ad- 
vantage of having a self-consistent treatment of the Compton 
up-scattering process. 

Several other models describe the process of Compton up- 
scattering of seed photons, with different degrees of complica- 
tion and various limitations. The models COMPTT and COMPST 
(see Titarchuk 1994, and references therein) describe a thermal 
Comptonization in a static plasma so, despite the possibility of 
obtaining statistically satisfactory results in terms of X-ray spec- 
tral fitting, they are not suitable to be applied to accreting sources 
at highly sub-Eddington rates as RXJ0440. 9+4431 (in these 
cases the accreting material is expected to be nearly in free-fall). 

The model COMPTB (Farinelli et al. 2008) has been pro- 
posed to describe the emission of low-mass X-ray binaries 
for both the cases of simply thermal or thermal plus bulk 
Comptonization. However, the correlation between the electron 
temperature kT e , the spectral index a and the bulk parameter S 
does not yield an independent determination of the whole pa- 
rameter set. Assuming that bulk Comptonization is an impor- 
tant effect, the model implementation requires that the param- 
eter kT e is fixed or allowed to vary in a very narrow range, 
inferred from other considerations. It is also important to note 
that the COMPTT, COMPPS, and COMPTB models do not consider 
the presence of a magnetic field so that the classical Thomson 
or Klein-Nishina electron scattering cross-sections are assumed. 
The Be X-ray binary nature of RX J0440. 9+4431, together with 
its relatively high pulsed fraction (Usui et al. 2012), suggests 
that the NS in this system is endowed with a magnetic field 
significantly higher than in the low-mass systems. This is ex- 
pected to funnel efficiently the accreting material onto the neu- 
tron star's polar caps, leading to an electron temperature that is 
hardly predictable a priori. 

The model proposed by Becker & Wolff (2007) is also 
not appropriate to describe the X-ray emission observed from 
RXJ0440. 9+4431, as it assumes a velocity profile for the in- 
flowing plasma in which an intense radiation field originates 
a collision-less shock above the NS surface. In this case, 
the radiation creeps through the lateral walls of the accretion 
stream giving rise to the so-called “fan beam”. In the case of 
RX J0440. 9+443 1 , the emission pattern from the NS is expected 
to assume rather the shape of a “pencil beam”, as at low lumi- 
nosity the radiation is believed to be able to escape directly in 
the direction of the accretion stream and no dramatic changes in 
the pulse profile could be revealed over two orders of magnitude 
in luminosity (La Palombara et al. 2012; Usui et al. 2012). The 
latter assumption can be roughly demonstrated in the present 
case by comparing the magnitude of the gravitational force of 
the neutron star and the radiation pressure from the polar cap 
where matter is thermalized. We find that 

Jfgrav ~ 1.3 x 10 14 mr 6 _I , (5) 

where m and r& are the neutron mass in solar masses and the 
distance from the surface in units of 10 6 cm, respectively. On the 
other hand, the radiative force per unit mass is given by 

,^ rad * 1 .4 x H)' 3 kT m A f(B), (6) 

where kT m is the temperature of the thermal mound and 
f(B ) < 1 is a correction factor to the Thomson cross-section 
due to the presence of the magnetic field. For photons propa- 
gating along the vertical direction, thus almost parallel to the 
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Table 4. Best fit parameters of the COMPMAG model using a free-fall 
velocity profile of the accreting matter. 



April 2010“ 

September 20 10* 

A h [ 10 22 cm -2 ] 

0.36 + 0.01 

0.28 ± 0.02 

kT B B [keV] 

1.50 + 0.01 

1.32 + 0.01 

kT e [keVl 

13+ 1 

62+1 

T 

0.12+ 0.01 

0.09+ 0.01 

A> 

>0.56 

0.44 ± 0.01 

7*0 

0 30 +0 05 

o.JW-002 

0 32+ 0 01 
-0.04 

N 

112 + 5 

388 + I 

^Gauss [keV] 

6.33 ± 0.03 

- 

A Gauss 

(6.59 + 1.2) x 10“ 4 

- 

FluX 0 .l_ 200 keV C 

3.4 

3.0 


1.21/569 

0.93/558 


Notes. (a) Observation A in Table 2. (6) Observation S in Table 2. 
(c) Estimate of de-absorbed flux in units of 10 -9 ergcm -2 s _1 , referred 
to Swift/XKT observation. 


magnetic field, the Thomson cross-section is strongly reduced, 
with f(B ) <k 1 (Becker & Wolff 2007). For a typical neutron star 
mass and using the best fit values of kT m reported in Tables 3 
and 4, it is evident that ^ rad <s 

The model of Becker & Wolff (2007) was extended recently 
(COMPMAG; Farinelli et al. 2012). In this new implementation, 
it is possible to tune the velocity profile of the in-falling material. 
A broad-band spectral coverage is however necessary to obtain 
suitable constraints in the physical parameters of COMPMAG, and 
thus the single PCA spectra could not be used for this purpose. 
More satisfactory results were obtained in the multi-instrument 
broad band spectral analysis of the outbursts in April 2010 and 
September 2010. We have chosen to use a free-fall like profile 
with fixed parameter rj = 0.5, we have also fixed the albedo 
at the NS surface to one, while the other parameters are de- 
termined by the spectral fit. They are the flow velocity at the 
NS surface (j3o, expressed in unities of c), the accretion col- 
umn radius (ro, expressed in Schwarzschild radii), the seed black 
body temperature {kT% b), the electron temperature in the accre- 
tion stream ( kT e ), their opticsl depth (r), and the model normal- 
ization ( N ). A summary of these results is reported in Table 4. 
The temperature of the seed thermal photons is in fair agreement 
with the BMC model. The increase of the electron temperature 
at decreasing luminosity can be identified as the origin of the in- 
creased cut-off energy, while the normalization, N, of COMPMAG 
is defined such that m = N 1,2 D io, and thus gives an estimate 
of the apparent blackbody emitting area. It is worth noting that 
COMPMAG assumes a uniform blackbody photon field distribu- 
tion across the accreting column, so the inferred value R ^ does 
not actually represent the hot spot dimension on the neutron star 
surface, but rather the total area obtained by emission at differ- 
ent heights above the neutron star surface (so that in general 
Rkm > ro). This is confirmed here: we obtain Rkm ~ 6.5 and 
Rkm ~ 3.5 for the first and second observation, respectively (at 
a source distance of 3.3 kpc), while the column radius remains 
constant within the uncertainties at about 1 km. 

We note that the absorbing column as measured by 
XMM-Newton / PN is significantly higher than the Swift/XKT de- 
terminations, while the spectral parameters of the BMC slightly 
deviate from the trend derived from the PCA data set (Fig. 4). 
This can be due to the different energy ranges over which spectra 


are measured which can lead to a difference in the determina- 
tion of the spectral slope and thus of the absorbing column. 
Furthermore, we note that according to Romano et al. (2005), the 
PN detector yields spectral parameter for the calibration source 
PSRB0540-69 which deviate from EPIC/MOS and Swift/XKT 
in the same direction as our data, i.e., higher Ah, although by 
a smaller amount. Therefore, we cannot exclude that the appar- 
ent discrepancy we have found is influenced by the calibration 
systematics. We have verified that using a different absorption 
model (tbabs) with the abundances of Wilms et al. (2000) does 
not remove the discrepancy and yields an equally acceptable fit 
with compatible parameter value, except the absorbing column, 
which is systematically higher than in the model used in this 
paper. We remark that the absolute value of the absorbing col- 
umn is of marginal importance for the results discussed here, 
and that the relatively low-resolution spectra analyzed here do 
not allow us to distinguish between the models. The values we 
have derived are similar with the one determined by Tsygankov 
et al. (2012), while Usui et al. (2012) use just RXTE data (above 
3 keV) and get absorbing column roughly one order of magni- 
tude in excess: we have verified that excluding the Swift/XKT 
data, we would have obtained similar values, which we conclude 
to be inaccurate as they are hampered by the lack of low-energy 
coverage. 

Among all these different spectral models, the BMC is able 
to provide a self-consistent and robust physical description of 
the accretion process in RXJ0440. 9+4431 over the different 
luminosity levels even exploiting a limited energy band. The 
COMPMAG model is able to provide complementary informa- 
tion on the accretion process, which are not present in the 
BMC model, although with more structured assumptions on the 
geometrical properties of the accretion stream, in particular 
a significant vertical extent of the emitting plasma above the 
NS surface. Flowever, it generally requires a large energy cover- 
age and it is of limited use when data are not available to provide 
constraints on the source emission above ~20keV. 

4.3. Spectral parameters 

The observations of RX J0440. 9+443 1 in 2010 and 2011 al- 
lowed us to study spectral changes in the X-ray emission from 
the source over two orders of magnitude in X-ray flux and mass 
accretion rates (~10 14-16 g s -1 ). To compare the results of our 
spectral analysis at different accretion rates, we plot in Fig. 6 a 
selection of our BMC best-fit unfolded spectral models obtained 
at increasing flux levels from the fits to the RXTE data. We note 
that the spectra of the lower persistent state are rather soft, they 
become harder at intermediate states, and display a more pro- 
nounced curvature at the highest luminosities. This is in agree- 
ment with the measurement of a high-energy exponential cut-off 
at ~30keV in the RXTE/PCA data during the two observations 
performed when the source flux was >10 -9 ergcm _2 s _1 . The 
photon index of the Comptonization Green’s function slightly 
decreases in comparison to the observations performed at the in- 
termediate flux level, also indicating a hardening of the source 
X-ray emission. At the intermediate state, INTEGRAL/ISGRI 
data were available to provide coverage at the higher energies, 
an exponential cut-off was measured at (53 4 [g)keV. This sug- 
gests that the cut-off energy increases at lower luminosities so 
that the hard X-ray spectrum presents the hard tail expected 
when bulk Comptonization is the dominant reprocessing chan- 
nel (Becker & Wolff 2005). Observations performed at the lower 
fluxes (<10 _1 ° ergcm _2 s _1 ) showed a steeper Comptonization 
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Fig. 6. Selection of unfolded RXTE/PCA spectra at different luminosity 
states based on the BMC model. The corresponding observation ID is 
also indicated. 


index and a cooler thermal photon temperature, they are thus 
characterized by a significantly softer emission. These results 
are confirmed by the model-independent analysis of the source 
emission provided by the HID in Fig. 5. The observations in 
which a cut-off energy of 30 keV was measured are represented 
with the two blue diamonds located in the upper right corner of 
the plot. The PCA observations performed at intermediate levels, 
for which quasi-simultaneous INTEGRAL/ISGRI data were also 
available, are represented with the two green stars located im- 
mediately below the blue diamonds. The source emission hard- 
ens remarkabley at count rates ~80-100cts/s, i.e., where the 
change in the cut-off energy is also revealed by the spectral fit- 
ting. The flux corresponding to the lowering of the cut-off energy 
is 1 .3 + 1 .8 x 10 -9 erg cm -2 s -1 , and translates into a luminosity 
~2 x 10 36 erg s -1 . 

According to the discussion of Becker et al. (2012), at low 
luminosity ( L x < 10 34-35 erg s _1 ) matter strikes the NS surface 
after passing through a gas-mediated shock. At higher luminos- 
ity, the radiation pressure begins to affect the accretion process 
so that a radiative shock forms, and plasma deceleration takes 
place through Coulomb interactions in an vertically extended at- 
mosphere. At this stage, the efficiency of the thermal cooling 
of the high-energy electrons responsible for the Comptonized 
emission increases significantly and can produce the observed 
change in the cut-off energy. However, the luminosity level at 
which this transition happens is not well predicted on theoretical 
basis and is here tentatively identified on observational grounds 
at ~2 x 10 36 ergs _1 . It is worth noticing that in all observa- 
tions, RX J0440.9+443 1 remains relatively faint in the X-rays 
compared to the giant outbursts of transient Be/X-ray binaries 
(Lx ~ 10 38 erg s -1 ). As a consequence, the HID does not show 
the typical upper branch turning towards softer colors typically 
observed during the brightest outbursts of these sources, which 
might be interpreted as the onset of the fully radiation dominated 
accretion column (Reig & Nespoli 2013; Becker et al. 2012). 

To test further our understanding of the accretion processes 
in RXJ0440. 9+4431 given above, we have estimated from the 
BMC best-fit parameters the terminal velocity of the accretion 
flow at the NS surface for all the RXTE spectra and reported our 
results in Fig. 6. For this purpose, we follow Eqs. (91), (92) of 
Becker & Wolff (2007) adapted to the notation used in the previ- 
ous sections of our paper. By equating the energy flux dissipated 



Fig. 7. Luminosity dependency of the velocity of the accreted matter at 
the NS surface obtained from Eq. (10) assuming C = 0.25. 


by the thermally emitting surface to the flux of kinetic energy, 
we obtain 

^sbL 4 b = l -Jv 2 , (7) 

where cr SB is the Stefan-Boltzmann constant, L BB the black body 
temperature, v the stream velocity, and 



( 8 ) 


The radius of the hot spot emitting the thermal radiation, r BB , can 
be derived from the measured spectral parameters using Eq. (2). 
We use also the standard conversion between mass accretion rate 
and X-ray luminosity: 


M = 


4 nD 2 F x R 2 t 
GM, 


(9) 


where D is the source distance, F x is approximated by L2-30kev, 
G is the gravitation constant, R, and M, are the NS radius and 
mass. By taking into account that Abmc = Lig/D 2 0 (see Eq. (2)) 
and expressing the free-fall velocity at the NS surface as vs = 
V2 GM,/R„ we obtain that 


v_ _ I Nbmc 
m ~ V 4.8 x 10 7 CF X ’ 


(10) 


where C is the geometrical factor of Eq. (2) and F x is expressed 
in ergcm _2 s _1 . By assuming that the thermal emitting surface 
is a circular slab (C = 0.25), we derive a terminal velocity 
of the accretion stream which is ~0.85uff for all the observa- 
tions, apart for the two corresponding to the highest luminosities 
(Fig. 7). In these cases the terminal velocity marginally exceeds 
the free-fall limit. This might indicate that the thermal mound is 
building-up on the NS surface so that the geometrical factor C 
of the thermally emitting region increases and approaches that 
of the spherical geometry case (C = 1). This indication agrees 
with the interpretation that the spectral change at the peak of 
the flaring activity of RXJ0440. 9+4431 is due to the formation 
of a radiatively dominated accretion column suggested by the 
decreasing cut-off energy. We note that the above calculations 
agree fairly well with the base velocity of the COMPMAG best-fit 
parameters (flo ~ 0.5 ~ 0.8uff/c). 

Throughout the whole luminosity range, the variation 
of the parameters of the thermal emision component in 
RXJ0440. 9+4431 could be satisfactorily described by the re- 
lations: r BB cc F’0-39±o.02 ; kTm K f o.070±o.oos ( see Sect 3 3 and 
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Fig. 4). In the case of accretion onto a neutron star from the stel- 
lar wind, it is expected that the radius of the hot spot is inversely 
correlated to the luminosity, as at higher accretion rate matter 
should penetrate along magnetic field lines which are more and 
more concentrated towards the NS magnetic poles (Arons & Lea 
1980). Conversely, if the transfer of matter is mediated by an ac- 
cretion disk, the area of the hot spot is proportional to the ra- 
tio /^Aifvenoc if 1 ( see > e 'S-’ Lamb et al. 1973; White et al. 
1983). The correlation determined from the fits to the data thus 
support the idea that RX J0440. 9+443 1 is continuously accreting 
from a surrounding disk (despite the opposite would be expected 
for the persistent activity of a Be/X-ray binary, see Sect. 4.1). 
However, the measured relation between r B B and the source lu- 
minosity is much steeper than expected in the case of disk accre- 
tion 8 * . This might indicate either that matter penetrates through 
closed magnetic field lines at the border of the magnetosphere 
or that the neutron star magnetic field structure is a more com- 
plicated than a simple dipole close to the surface. The former 
idea was suggested to justify the magnetospheric instabilities ob- 
served in the outbursts of EXO 2030+375 (Klochkov et al. 2011) 
and A 0535+262 (Postnov et al. 2008). The latter has been in- 
voked, e.g., by Ferrigno et al. (2009) to interpret the soft com- 
ponent in the broad-band spectrum of 4U01 15+63 and can be 
related also to accretion-induced deformations of the magnetic 
field and instabilities in the accretion mound which are shown 
to exist on the base of recent numerical simulations (Mukherjee 
et al. 2013). 

The way in which the accreting matter penetrates the 
NS magnetosphere and accretes onto the compact star surface is 
also producing significant changes in its spin period. In Sect. 3.2, 
we have shown that the spin period values measured from the 
observations of RXJ0440. 9+4431 are not precise enough to 
constrain the source orbital properties and/or evaluate accretion 
torques. However, we note that RXJ0440. 9+4431 showed a re- 
markable spin-down over the past 12 years from 202 s (Reig & 
Roche 1999) to ~206s. If we assume that RX J0440. 9+443 1 is 
accreting from a disk at a non negligible rate, then the measured 
spin down can be interpreted in terms of coupling between the 
neutron star magnetic field lines and the external regions of the 
accretion disk located beyond the so-called corotational radius 
(Wang 1995, 1996). This assumption has been sometimes used 
to constrain the NS magnetic field. As pointed out by Bozzo 
et al. (2009), the uncertainties on the location and the extension 
of the boundary layer in the disk accretion models for magne- 
tized neutron stars are relatively large. Therefore, we could just 
derive that the magnetic field for the NS in RX J0440. 9+443 1 is 
of the order of few x 10 12 G. This is in agreement with the timing 
signatures of the broad-band noise discussed in Tsygankov et al. 
(2012). The non-detection of the previously reported cyclotron 
scattering absorption line discussed in Sect. 3.3 is not a counter 
argument to this magnetic field estimate, as these features are not 
ubiquitous in the spectra of accreting magnetized neutron stars. 


5. Conclusions 

In this work, we reported on the analysis of the last outbursts 
detected from the persistent Be X-ray binary RXJ0440. 9+4431, 
and used all archival observations of this source performed in 
2010 and 2011 to investigate its spectral properties in a broad 
X-ray energy domain (0.3-150keV). 


8 To infer the thermal component radius from the model parameters, 

we have assumed a constant geometry (parameter C of Eq. (2)). 


We obtained the most accurate measurement of the source 
orbital period so far at (150.0 ± 0.2) days and performed a de- 
tailed analysis of the source spectrum covering more than two 
orders of magnitude in luminosity. We used the bulk motion 
Comptonization model (BMC, Titarchuk et al. 1997), which pro- 
vided a self-consistent picture of the accretion processes in 
RXJ0440. 9+4431. In particular, we found a coherent depen- 
dency of the size and temperature of the thermal emitting re- 
gion as a function of the luminosity. These quantities do not 
follow the relations expected for accretion onto strongly mag- 
netized neutron stars, and we argue that either accreting mat- 
ter penetrates through closed magnetic field lines at the bor- 
der of the magnetosphere or that the neutron star magnetic field 
configuration deviates from a dipole close to the surface. 

Analyzing the spectral changes at the hard X-ray band 
(>10 keV), we have identified the luminosity level at which 
the accretion flow dynamics begin to be affected by radiation 
pressure (L x ~ 2 x 10 36 ergs _1 ). We estimated the in-falling 
plasma velocity close to the neutron star, which is consistent 
with the predictions from a physically motivated model. This 
provided an independent indication of the building-up of a ra- 
diative accretion column at the highest luminosity reached by 
RXJ0440. 9+4431 from the likely bending of the thermal emis- 
sion surface. RXJ0440. 9+4431 was particularly well suited 
for this study because the variations of its thermal emission 
could be unambiguously revealed during the higher and lower 
X-ray active states. This is not the case for other Be X-ray bi- 
naries (see, e.g., the controversy on the spectral modeling of 
EXO 2020+375, Reig & Roche 1999; Klochkov et al. 2007; 
Wilson et al. 2008). 

We compared our findings on the accretion process in 
RXJ0440. 9+4431 by exploring the most recent spectral mod- 
els accounting for thermal emission and its Comptonization in 
a strong magnetic field. The COMPMAG model could provide a 
number of further details for the interaction between the neutron 
star hosted in this source and the in-falling plasma, but the us- 
ability of this model is limited by the statistics of the available 
data. More sensitive missions will hopefully cover the required 
broad energy range (2-50keV) (e.g., LOFT, Feroci et al. 2012) 
and will allow us to study in greater detail the accretion process 
at the lower end of the luminosity range of Be/X-ray binaries by 
means of complex physically motivated models, which cannot 
be significantly constrained at present. 
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